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Abstract 

The  sound  field  generated  by  an  open  rotor  and  that  radiated  by  the  same  rotor  placed  inside  a semi- 
infinite duct  is  simulated  using  a Computational  Aeroacoustic  (CAA)  multi-domain  method.  The 
linearized  Euler  equations  in  cylindrical  coordinates  are  used  as  governing  equations  in  solving  this 
problem.  Both  Tam's  radiation  and  PML  boundary  conditions  are  used.  Comparisons  of  results  by  using 
the  two  different  far-field  boundary  conditions  will  be  given.  A multi-domain  Cartesian  grid  system  is 
used  so  that  the  number  of  grid  points  can  be  kept  as  a minimum.  A combination  of  DRP  and  a specially 
optimized  cell-centered  high  order  differencing  scheme  is  implemented  in  the  area  of  two  Cartesian  grid 
interfaces.  The  effect  of  mean  flow  including  a shear  layer  on  the  noise  radiation  off  the  ducted  rotor  is 
studied. 

1.  Introduction 

The  prediction  of  open  or  ducted  rotor  noise  as  shown  in  Fig.fl  ) is  a challenging  problem  as  there  are 
large  difference  in  scales  between  the  noise  source  near  the  rotor  and  the  sound  field.  The  most 
commonly  used  methods  for  this  type  of  problem  are  integral  methods  based  on  both  FW-H  and 
Kirchhoff  formulations[l,2].  These  approaches  to  compute  noise  radiation  from  rotors  in  the  subsonic 
flow  regime  are  split  into  two  steps.  In  a first  step,  the  near-field  aerodynamics  and  acoustics  close  to  the 
source,  or  on  the  surface  when  using  FW-H  with  impenetrable  surface,  are  computed  using  an 
aerodynamics  code.  This  then  serves  as  input  to  the  integral  method  which  computes  the  acoustic 
pressure  at  any  desired  observer  location  in  the  far-field.  For  the  computation  of  acoustic  propagation 
over  large  distances  without  dissipation,  these  hybrid  methods  arc  much  more  computationally  efficient 
than  CFD  methods,  but  they  suffer  from  the  limitation  that  acoustic  waves  are  assumed  to  propagate 
through  a medium  which  is  largely  at  rest  or  at  most  in  uniform  motion  relatively  to  the  observer. 

In  the  present  paper,  the  sound  field  generated  by  an  open  and  ducted  rotor  is  simulated  using 
Computational  Aeroacoustics  (CAA).  The  CAA  method  [3]  has  the  advantage  over  a CFD  method  in 
accurately  predicting  small  amplitude  acoustic  fluctuations  and  their  correct  propagation  to  the  far  field 
even  in  non-uniform  flow.  The  linearized  Euler  Equations  (LEE)  are  used  as  governing  equations  in 
solving  this  problem.  The  problems  considered  are  taken  from  the  third  CAA  benchmark  problems 
workshop[4].  In  this  benchmark  problem,  a rotor  is  simplified  as  distribution  of  a rotating  force  acting  on 
the  fluid.  In  the  far-field  area  of  the  computation  domain  as  well  as  the  duct  inflow  in  the  ducted  rotor 
case,  both  Tam's  radiation  and  Perfect  Matched  Layer  (PML)[5]  boundary  conditions  are  used. 
Comparisons  of  results  by  using  the  two  different  far-field  boundary  conditions  will  be  given.  Because  of 
the  disparity  between  the  length  scales  in  the  source  region  and  the  acoustic  wavelengths,  a multi-domain 
Cartesian  grid  system  is  used  in  which  local  resolution  enhancement  is  achieved  by  stepping  the  grid  size 
of  neighboring  Cartesian  grid  blocks  by  factors  of  two  so  that  the  number  of  grid  points  can  be  kept  as  a 
minimum.  In  the  area  of  two  Cartesian  grid  interfaces,  a combination  of  Tam  & Webb’s  DRP  and  a 
specially  optimized  cell-centered  high  order  differencing  scheme  is  implemented  [6,7], 

The  report  is  structured  as  follows.  The  governing  equations,  numerical  schemes,  boundary  conditions 
and  multi-domain  method  are  discussed  at  first.  The  noise  radiation  from  open  and  ducted  rotor  with 
subsonic  and  supersonic  rotational  speeds  are  studied.  Some  of  the  results  are  compared  with  the 


Paper  presented  at  the  RTO  AVT  Symposium  on  "Ageing  Mechanisms  and  Control: 
Part  A - Developments  in  Computational  Aero-  and  Hydro- Acoustics  ” , 
held  in  Manchester,  UK,  8-11  October  2001,  and  published  in  RTO-MP-079(1). 
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analytical  solutions.  In  order  to  demonstrate  the  effect  of  mean  flow  on  the  noise  radiation  off  the  ducted 
rotor,  a mean  flow  is  introduced  including  a shear  layer,  separating  the  ambient  medium  from  the  jet 
flow  exiting  out  of  the  duct. 

2.  Governing  Equations 

The  basic  set  of  equations  for  the  present  analysis  are  the  three  dimensional  LEE  with  a non-uniform 
mean  flow  in  cylindrical  coordinates.  Under  the  assumption  that  the  perturbations  are  small,  LEE  can  be 
obtained  by  splitting  the  flow  variables  of  the  Euler  equations  into  the  mean  and  perturbation  terms.  The 
linear  perturbation  field  represents  acoustic,  vorticity  and  entropy  waves.  The  linearized  perturbation 
equations  or  LEE  about  the  mean  flow  in  cylindrical  coordinates  then  have  the  following  forms  : 
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In  which  u , v'  and  w'  are  the  axial  (x),  radial  (r)  and  azimuth  ( 0 ) velocity  perturbation  components 
respectively,  p'  and  p'  are  the  density  and  the  pressure  respectively.  H,v,w,p  and  pare  the 
corresponding  mean  values,  and  the  subscripts  x,  r and  0 denote  their  partial  derivatives.  The 
( Fx , Fr , Fg ) denote  the  components  of  the  source  term.  All  quantities  are  written  in  non- 
dimensionalized  form  with  respect  to  length  scale  b (radius  of  the  blade),  velocity  scale  c (sound  speed), 
timescale  b /c,  density  scale  pn  (ambient  density) , and  pressure  scale  p0c2. 


When  a parallel  mean  flow  condition  (m*0,v  = w = 0 and  p = yp  = 1 ) is  introduced,  the  Eq.(l)  can  be 
simplified  as 
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Since  for  the  parallel  mean  flow  there  is  no  coupling  of  the  velocity  and  pressure  perturbation  with  the 
density,  the  density  equation  contained  in  Eq.(l)  is  neglected  in  Eq.(2). 

To  simplify  the  problems  while  keeping  the  main  physics,  the  rotor  is  replaced  by  a distribution  of 
rotating  body  forces.  The  components  of  body  force  or  the  source  terms  (Fx,Fr  ,Fg  ) in  Eq.(l)  and 
Eq.(2)  are  given  from  the  third  CAA  benchmark  problems  workshop  as  follows. 


Fx(x,  r,9,t) 

Fx(x,r) 

Fr(x,r,6,t) 

— Re- 

0 

Fe(x,r,6,t) 

Fg(x,r) 

where'  Re{  }is  the  real  part  and  (Fx,Fe ) denotes  the  body  force  distribution  in  r and  x.  In  addition,  the 
m is  the  number  of  rotor  blades  and  Q is  the  angular  velocity  of  the  rotor.  Therefore,  the  frequency  of 
the  radiated  sound  is  mQ . 


The  actual  rotor  noise  is  a 3-D  problem,  but  by  factoring  out  the  azimuthal  dependence  with 
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the  problem  can  be  reduced  to  a 2-D  problem.  Then  the  governing  equations  for  (n,v,w,/5)  can  be 
obtained  by  inserting  Eq.(3)  and  Eq.(4)  into  Eq.(2), 
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The  Eq.(5)  will  be  used  to  solve  open  and  ducted  rotor  problems.  The  solutions  to  the  Eq.(5)  can  be 
directly  discretized  in  x and  r for  Cartesian  grids  and  a given  parallel  mean  flow.  The  DRP  finite 
difference  scheme[8]  and  the  artificial  selective  damping  [9]  are  used  for  the  spatial  discretization. 

2.  Boundary  Conditions 


Due  to  the  limitation  of  a finite  computational  domain,  numerical  outer  boundary  conditions  are  needed 
for  the  acoustic  simulation.  At  the  outer  boundary  of  the  computational  domain,  two  types  of  non- 
reflecting boundary  conditions  are  implemented.  These  are  asymptotic  radiation  and  outflow  boundary 
conditions  as  well  as  PML[5]  conditions.  The  details  of  the  radiation  and  outflow  boundary  conditions 
developed  by  Tam  and  Dong  were  given  in  ref.[10].  In  the  case  of  ducted  rotor,  the  solid  wall  boundary 
condition  must  be  satisfied  at  the  surface  of  the  duct  wall.  In  addition,  along  the  axis  of  the  cylindrical 
coordinates;  i.e.,  r = 0 , a special  set  of  equations  are  used. 
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2.1  PML  Boundary  Conditions 


The  asymptotic  radiation  and  outflow  boundary  conditions  developed  in  ref.  [10]  can  only  be  used  for 
weakly  non-uniform  flow.  Its  use  for  the  flow  with  strong  shear  layer  across  the  far-field  boundary  area 
is  not  suitable.  PML  conditions  arc  implemented  as  an  absorbing  boundary  condition.  One  advantage  of 
using  PML  is  that  PML  can  be  easily  implemented  for  non-uniform  mean  flow.  The  PML  technique  has 
proved  to  be  effective  and  easy  to  implement.  For  the  open  and  ducted  rotor  problems,  the  PML 
equations  used  in  the  PML  domain  are  derived  by  splitting  the  variables  of  the  governing  equations  (5)  : 
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where  u - w,  + u2 , v = v,  + v2 , w = w,  + w2  and  p = pt+p2-  The  absorption  coefficients  <JX  and  a r 
are  varied  gradually  in  the  PML  domain.  As  shown  in  Fig.f3a~)  and  Fig.(8a')  the  PML  domains  used  for 
the  open  and  the  ducted  rotor  case,  Ox  is  set  zero  at  the  right  and  left  PML  domains  and  ar  is  zero  at  the 
top  PML  domain.  At  the  comer  regions,  both  coefficients  are  nonzero.  In  the  present  paper,  the 
following  form  has  been  used  for  the  variation  of  <TV  or  <J r within  the  PML  domains: 

CT  = (Tmax(1.0-e-(‘,/D)2)  (7) 

in  which  D is  the  width  of  the  PML  domain,  d is  the  distance  from  its  interface  with  the  interior  domain 
and  <7max  is  the  maximum  value  of  G , which  is  dependent  on  grid  size. 


2.2  Solid  Wall  Boundary  Condition 

For  the  ducted  rotor,  on  the  surface  of  the  duct  wall,  the  fundamental  physical  requirement  is  the  non- 
penetration boundary  condition: 

v»/7=0  or  unx  + v«^  + wne  =0  (8) 

where  n is  the  unit  normal  on  the  wall  and  can  be  written  as  h = (n f^n0)  = (0,  1,  0)for  a plane 
wall.  Then  Eq.(8)  becomes 


v = 0 (9) 

By  taking  the  dot  product  of  the  momentum  equations  of  Eq.(5)  with  n for  the  plane  wall,  the  Eq.(8)  is 
equivalent  to  : 


Sp_ 

dt 


(10) 


(SYA)  2-5 


In  the  case  of  the  ducted  rotor,  both  Eq.(9)  and  (10)  are  used  on  the  wall  surface.  For  the  flat  wall,  Eq.(9) 
can  be  easily  implemented.  To  satisfy  Eq.(lO),  the  numerical  treatment  in  conjunction  with  high-order 
finite  difference  schemes  derived  by  Tam  and  Dong  [9]  is  implemented.  In  this  method,  ghost  values  of 
pressure  are  chosen  so  that  Eq.(10)  can  be  satisfied  on  each  boundary  condition  enforcement  point. 

2.2  Conditions  at  the  Centerline  r-  0 

As  r — » 0 , Eq.(5)  has  a numerical  singularity  and  cannot  be  used  directly.  Since  for  m & 0 on  the 
centerline  r — 0 , v = w = p = 0 , the  singularity  terms  in  Eq.(5)  may  be  replaced  by  partial  derivative 
terms  using  L’Hospital’s  Rule,  which  are 


at  r = 0 . 


v _ 3v  w _ dw  p _dp 
r dr  ’ r dr  * r dr 


(11) 


3.  Multi-Domain  Grid  Design 

Because  of  the  disparity  between  the  length  scales  in  the  source  region  and  the  acoustic  wavelengths,  a 
multi-domain  Cartesian  grid  system  is  used.  The  computation  domain  is  divided  into  several  sub- 
domains  enclosed  by  rectangular  solid  lines  as  shown  in  Fig.t2al  for  an  open  rotor  computation.  For 
clarity,  only  part  of  the  outside  domain  is  given  in  Fig,f2al.  The  highest  resolution  is  used  in  the  area 
closest  to  the  rotor  blade  or  the  source  region  respectively.  When  moving  away  from  source  region  to  the 
acoustic  field,  the  grid  size  is  stepped  up  by  factors  of  two  from  inner  to  outer  sub-domains.  In  principle, 
an  overlapping  grid  technique  for  the  use  of  high  resolution  schemes  [11]  can  be  used  in  the  interface 
area.  To  make  the  computation  more  efficient,  a special  interface  treatment  between  Cartesian  grids  was 
used.  Since  the  grid  spacing  in  adjacent  Cartesian  areas  differs  by  a factor  of  2,  it  is  not  difficult  to  form  a 
7 points  DRP  stencil  at  every  point  in  region  A with  mesh  size  Ax  or  Arand  every  point  in  region  B 
(mesh  size  Ax/2  or  Ar/2  ) except  for  points  in  the  interface  region,  which  is  within  a distance  of  three 
points  from  the  interface  in  region  B as  shown  in  Fig.f2bl.  In  the  interface  region,  when  grid  lines  on  two 
grids  coincide  with  each  other,  a cell-centered  symmetric  stencil  with  optimized  coefficients  is  used,  such 
as  the  point  marked  'cc'  in  FigJ2b).  The  stencil  used  at  'cc'  can  be  written  in  the  form 


The  optimized  stencil  coefficients  are 

a_  3 = -a 3 = -0.0046875 
a_2=-a2  = 0.065 1 04 1 666667 
a,  = -a,  =-1.171875,  a0  = 0.0 

For  the  points  on  the  non-coincidence  lines,  function  values  on  several  ‘hanging’  points  marked  as  ‘h’  in 
Fig.(2bl  are  required.  The  values  on  these  ‘hanging’  points  are  explicitly  obtained  by  a two  dimensional 
Lagrange  interpolation.  The  interface  treatments  used  here  have  been  validated  in  ref.[6,7j  on  the  basis  of 
CAA  Benchmark  problems. 

4.  Artificial  Selective  Damping  (ASD) 

The  ASD  terms  are  added  in  the  right  hand  side  of  the  discretized  finite  difference  equations  in  order  to 
eliminate  spurious  short  waves  and  improve  numerical  stability.  In  the  interface  regions  among  Cartesian 
grids,  the  irregularities  are  created  by  using  different  grid  size  and  stencil.  These  irregularities  are  also 
very  efficient  in  generating  short  spurious  numerical  waves.  For  the  field  points  where  the  cell-centered 
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symmetric  stencil  Eq.(12)  is  used,  a damping  stencil  with  the  same  configurations  was  derived.  The 
optimized  stencil  coefficients  of  the  damping  operator  are 


d_  3 = d}  = -0.01092915,  d_2  = d2  = 0.06004774 
d_s  =dt  =-0.2991186,  d0  = 0.5 


In  addition,  for  the  grid  points  located  at  the  centerline  r = 0 , the  imposition  of  stronger  ASD  is  required 
since  there  is  a change  of  governing  equations  which  may  result  in  the  generation  of  spurious  short 
waves. 


5.  Numerical  Simulations  and  Discussions 

Numerical  results  on  noise  radiation  from  the  open  and  ducted  rotor  will  be  discussed  in  the  following 
sections.  Because  the  variables  (u,  v,w,p)  used  in  the  governing  equations  are  complex,  the  results 
presented  are  taken  from  Eq.(4)  for  zero  azimuth  angle  0=0.  The  directivity,  D(<p) , of  the  radiated 
sound  is  defined  as 


D(<P)  = R2p'2(R,<t>,6,t)  (13) 


Where  is  the  time  average.  The  definition  of  R and  0 is  given  in  Fig.(l)  for  the  open  rotor. 

For  the  ducted  rotor,  R is  measured  from  the  open  end  of  the  duct. 


5.1  Open  Rotor 


The  sound  generated  by  an  open  rotor  is  simulated.  Following  body  force  distributions (Fx,Fe)  are  used 
to  simulate  the  impact  of  the  rotor  blades  on  the  surrounding  air[4]. 


Fe(x,r)  = 


Fx(x,r)  = 


\F(x)rJm(amNr)  r<  1 

[0  r > 1 

\F(x)Jm  (amNr)  r<  1 

[0  r > 1 


F(x)  = exp{-  (In  2 )(x  / /?„.  )2 },  hw  = 0.1 


(14) 


where  Jm(  ) is  the  m th-order  Bessel  function,  amN  is  the  Nth  root  of  J'm{amN  )=0  and  hw  is  the  half- 

2.7V 

width  of  the  Gaussion-like  forcing  function.  For  a B blade  rotor,  the  body  force  is  spaced  in  — radians 

B 

apart  and  spins  with  rotor  rotational  speed  £2 . The  m is  chosen  as  the  number  of  the  rotor  blades  B . The 
choice  of  the  Bessel  functions  in  Eq.(14)  has  no  other  significance  than  making  the  analytical  solution 
simple.  The  8-blade  rotor  ( 2?  = 8 ) is  involved  in  the  computation.  amN  is  chosen  as  a8 , = 9.64742  in 

which  N = 1 so  that  the  force  in  the  radial  direction  corresponds  to  a mono-peak  function.  In  ducted  rotor, 
N has  the  meaning  of  the  radial  mode  number. 

The  whole  computational  domain  is  xxr  = [- 17,  17]x  [o,  17].  The  problem  is  symmetric  about  the  x- 
axis.  Thus  only  the  solution  in  the  upper  half  x-r  plane  needs  to  be  computed.  The  symmetric  condition 
is  applied  for  (w,w,p)and  the  anti-symmetric  condition  is  applied  for  v. The  choice  of  the  grid  size 
depends  on  the  acoustic  wave  length  and  hw . In  general,  4 grid  points  are  necessary  to  resolve  the  half- 
width of  the  forcing  function  and  the  minimum  of  8 grid  points  per  wave  length  has  to  be  used  in  the  far- 
field. 
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Two  rotational  speeds  are  used,  one  with  subsonic  tip  speed  £2  = 0.85  and  the  other  with  supersonic  tip 
speed  £2  = 1.15.  The  grid  sizes  of  the  smallest  domain  are  Axs  = 0.025  for  £2  = 0.85  and  Ax,  = 0.02  for 
£2  = 1.1 5 respectively.  The  time  advancing  scheme  used  throughout  tire  paper  is  the  classical  4 stage 
Runge-Kutta  scheme  with  a time  step  At  = 0.005  if  there  is  no  other  specification.  In  the  open  rotor 
computation,  the  mean  flow  u(r)  = 0 is  assumed.  The  width  of  the  PML  domain  is  chosen  as  to  contain 
10  grid  points. 

Fig.  13  ) shows  a snapshot  of  the  acoustic  pressure  contours  and  the  details  of  the  multi-domain  grid 
system  for  the  open  rotor  with  subsonic  tip  speed.  The  rectangular  box  in  the  plots  is  the  interface 
between  two  Cartesian  domains.  The  results  show  most  of  the  acoustic  radiation  is  concentrated  in  the 
direction  around  the  plane  of  rotation.  The  results  in  zoom  scale  in  Fig.(3b)  also  demonstrate  that  the 
contour  lines  go  through  grid  interfaces  very  smoothly  without  any  noticeable  reflection.  Fig. 14)  gives 
the  directivity  at  R = 15.0  obtained  computationally  in  comparison  with  the  exact  solution  [4]  for 
£2  = 0.85.  There  is  good  agreement  between  the  numerical  simulation  and  the  analytical  solution 
although  the  latter  was  obtained  with  the  assumption  of  R — > °°  . To  compare  the  effect  of  different  far- 
field  boundary  conditions,  the  result  with  Tam’s  asymptotic  far-field  boundary  condition  is  also 
presented  in  the  figure.  Only  little  difference  was  found.  In  addition,  the  comparison  on  the  directivity 
for  different  rotational  speed  in  Fig.151  shows  that,  as  expected,  the  radiation  from  the  rotor  with  the 
supersonic  tip  speed  is  stronger. 


5.2  Ducted  Rotor 

For  the  ducted  rotor,  the  governing  equations  Eq.(5)  are  again  used  together  with  solid  wall  boundary 
conditions  on  both  sides  of  the  duct  wall.  In  general,  the  effect  of  the  duct  wall  will  force  the  acoustic 
disturbance  to  be  continuously  reflected  back  and  thereby  cause  cancellations  and  reinforcements. 
From  the  standard  modal  analysis  of  duct  acoustic  theory,  the  in-duct  pressure  for  a cylinder  duct  can 
be  expressed  as  superposition  of  various  rotational  modes, 

OO 

p{r,G,x,t)~  X HAmNJm  (KmNr)exp[int(e  - Clmtj  + iKxx\  (15) 

m N=0 


where  m and  N are  the  circumferential  and  the  radial  mode  numbers  respectively.  In  general,  KmN  is 
called  the  radial  wave  number  and  can  be  obtained  by  satisfying  the  wall  boundary  condition.  For  a 
hard  walled  duct,  KmN  is  the  N th  root  of  J'm  {KmN  )=  0 . £2„,  is  the  rotational  speed  of  the  m ’s  mode. 
When  the  body  force  felt  by  the  B-blade  rotor  is  steady  (as  in  the  present  paper),  the  rotational  speed 
of  all  modes  is  equal  to  £2  and  the  circumferential  mode  number  m can  only  be  chosen  as  m-fiB , 

( /t  = 1,2,3,...),  therefore  only  the  acoustic  signal  with  frequency  = l&p-  _ ^ ^ js  generated,  f is 

called  rotor  blade  passing  frequency' (BPF).  _ 

From  Eq.(15)  it  is  known  that  whether  an  acoustic  mode  can  be  propagated  along  the  duct  (x)  or  not 
depends  on  the  axial  wave  number  Kx . For  cylindrical  ducts  with  uniform  mean  flow  at  a Mach 
number  M -u  , Kt  is  obtained  by 


— KM  + yjK2  - K2mN(l-M^) 
1 -M2 


(16) 


where  K = pBi 2 is  the  non-dimensional  wave  number  of  the  acoustic  disturbance.  The  modes  with 
real  Kx  are  the  propagating  modes  and  called  cut-on  modes.  The  modes  with  non-zero  imaginary  part 
of  Kx  are  exponentially  decaying  in  amplitude  along  the  duct  and  therefore  called  cut-off  modes. 
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From  Eq.(16),  it  can  be  derived  that  the  acoustic  wave  with  BPF  and  its  high  harmonics  f can  be 
propagated  through  the  duct  only  if  the  rotational  speed  £2  satisfies, 

(17) 

m 

where  the  right  hand  side  of  Eq.(17)  is  called  cut-off  rotational  speed  £2  in  the  present  paper. 

For  the  ducted  rotor  with  the  same  body  force  function  as  Eq.(  14) , the  Eq.(15)  can  be  simplified  since 
the  given  body  force  function  can  only  excite  one  rotational  mode, 

p(r,0,x,t)  = Jm  (KmNr)exp[im(8  - £2 t)  + iKxx]  (18) 

where  KmN  = amN  and  m = B as  well  as  77  = 1 as  specified. 

5.2.1  Effect  of  Uniform  Mean  Flow  on  Noise  Radiation  from  Open  End  of  Duct 

In  the  ducted  rotor  case,  the  same  body  force  function  Eq.(14)  as  for  the  open  rotor  is  used.  An  infinitely 
thin  duct  wall  is  placed  at  r = 1.0  along  -4.1  <jc<8.0.  PML  conditions  are  implemented  on  the  left 
side  of  the  duct  to  simulate  an  infinite  length  of  the  duct.  Symmetry  about  the  x-axis  is  assumed.  The  size 
of  the  computational  domain  is  [- 4. 1,  16.5]x[0,  8]. 

For  the  case  without  considering  any  mean  flow  or  M = 0,  £2C  is  equal  to  1.206  according  to  Eq.(17) 
when  setting  m = B = 8 , N = 1 , and  KmN  = amN  = as , = 9.64742  as  in  the  open  rotor  case.  The  two 
rotational  speeds  £2  used  in  the  open  rotor  case  are  all  smaller  than  £2C  = 1 .206  so  that  all  the  acoustic 
disturbances  are  cut-off  or  evanescent.  This  can  be  seen  as  well  by  looking  at  the  time  history  outside  the 
duct.  Fig.f6t  gives  the  pressure  time  history  sensed  at  R = 6.0  and  <p  = 52°  for  £2  = 1.15 . R is  measured 
from  the  center  of  the  duct  end  located  at  (*,  r)  = (8.0,  0.0).  It  shows  that  the  intensity  of  sound 
radiation  decreases  very  quickly  after  the  initial  transient  state.  Thus  there  is  no  sound  radiation  out  of 
the  open  end  of  the  duct  even  for  slightly  supersonic  tip  Mach  number.  This  also  proves  that  the 
environment  around  an  acoustic  source  plays  an  important  role  in  noise  propagation. 

From  Eq.(17),  it  can  be  seen  that,  £2C  will  be  reduced  if  mean  flow  is  involved.  Therefore,  with  changing 
the  mean  flow  Mach  number,  the  condition  for  £2  to  either  represent  a propagational  or  non- 
propagational  signal  will  also  be  changed.  For  example,  for  the  case  M =0.5  inside  and  outside  of  the 
duct,  £2C  reduces  to  1 .044  so  that  the  noise  generated  by  the  rotor  with  £2  = 1.15  will  radiate  from  the 
open  end  of  the  duct.  The  pressure  time  history  for  the  ducted  rotor  with  £2  = 1.15  under  uniform  mean 
flow  M = 0.5  is  given  in  Fig,(7)  at  the  same  position  as  in  Fig.(6).  The  results  indicate  that  a time 
periodic  state  has  been  reached  after  the  initial  transient  state  due  to  sound  radiation  from  the  duct  end. 
Fig-f8)  shows  the  pressure  contour  at  / = 40  and  its  enlarged  plot  close  to  the  grid  interface.  The  results 
prove  again  the  numerical  treatment  is  successful.  It  should  be  mentioned  that  since  the  pressure 
disturbance  propagates  in  the  duct  on  a helical  path,  the  pressure  contour  given  in  Fig.(8)  is  only  a cut  in 
the  three  dimensional  spiral  pattern.  Due  to  the  mean  flow  convection  effect,  the  spinning  mode  will 
propagate  much  faster  in  downstream  direction  of  the  rotor  than  upstream  of  the  rotor.  The  difference  of 
the  length  scale  upstream  and  downstream  of  the  rotor,  as  shown  in  Fig.(8)  demonstrates  the  mean  flow 
convection  effect. 

5.2.2  Effect  of  Shear  Layer  on  Noise  Radiation 

In  previous  section,  it  was  demonstrated  that  the  effect  of  a mean  flow  in  the  duct  may  cause  the 
propagation  of  acoustic  modes  which  are  cut-off  for  vanishing  mean  flow.  Furthermore,  when  there  are 
gradients  in  the  mean  flow  typically  due  to  the  shear  layer  between  the  ambient  medium  and  the  jet  flow 
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emanating  from  the  duct,  instability  waves  may  be  generated  by  acoustic  excitation.  Finally,  the 
propagation  of  the  acoustic  waves  through  the  shear  layer  will  change  the  acoustic  directivity  pattern. 

In  order  to  show  the  main  effect  of  the  presence  of  such  a shear  layer  an  analytical  shear  flow  was 
introduced  separating  the  flow  field  in  r < 1 from  the  ambient  medium.  The  parallel  mean  velocity 
profile  is, 


u(r)  = 


M„  +(M 
M 


f(ta2)(T 


hfi  r>b 
h>r>0 


(16) 


where  the  mean  flow  is  uniform  inside  r < 1 with  iff  = 0.5  . The  Gaussian  distributed  function  is  used  to 
describe  the  variation  of  the  shear  flow.  The  half  width  of  this  function  is  hw  = 0.25 . The  other 
parameters  are  6 = 1.0  and  M.  = 0.0 . The  computation  domain  as  well  as  the  grid  arrangement  is  the 
same  as  in  Fig.(8).  PML  equations  of  Eq.(6)  are  again  used.  It  should  be  mentioned  that  the  shear  layer 
chosen  here  is  for  the  purpose  of  demonstrating  the  effect  of  the  shear  layer  on  noise  propagation  and  not 
to  show  the  excitation  of  the  instability  wave.  No  instability  mode  is  indeed  excited  under  the  considered 
signal  frequency. 

Fig.(9)  shows  the  instant  pressure  contour  at  two  different  times.  Apart  from  the  acoustic  waves  radiated 
from  open  duct  end,  depicted  in  Fig.f9al  a hydro-dynamical  instability  in  form  of  a wave  packet  is 
observed  inside  the  shear  layer.  The  instability  is  convected  out  of  the  computation  domain  at  later  times 
as  shown  in  Fig.f9b).  Therefore  these  instability  waves  are  believed  to  be  excited  purely  by  the  multi- 
frequency transients  occurring  when  the  force  Eq.(14)  is  turned  on  at  t = 0 . This  demonstrates  that  the 
considered  shear  flow  is  stable  for  the  acoustic  excitation  from  the  open  duct  end  at  £2  = 1.15.  The 
comparison  of  the  directivity  at  R = 6.0  for  uniform  and  shear  layer  mean  flow  is  given  in  Fig.(lO) 
which  shows  that  the  radiation  pattern  has  been  shifted  upstream  after  passing  through  the  shear  layer. 


6.  Summary 

The  noise  radiation  from  an  open  and  ducted  rotor  is  simulated  using  a CAA  multi-domain 
method.  The  results  for  the  ducted  rotor  show  that  there  is  no  sound  radiation  out  of  the  open 
end  of  the  duct  when  the  rotational  speed  of  the  ducted  rotor  is  below  the  cut-off  rotational 
speed.  The  sound  field  generated  by  an  open  and  ducted  rotor  can  be  extremely  different 
because  the  presence  of  the  duct  wall  will  cause  cancellations  and  reinforcements  of  both 
propagating  and  reflecting  waves.  The  effect  of  mean  flow  in  the  duct  may  cause  the 
propagation  of  the  acoustic  modes  which  are  cut-off  for  no  flow.  The  validations  using 
benchmark  problems  for  the  open  rotor  case  also  show  that  the  presented  method  has  good 
accuracy  in  terms  of  the  sehemes  and  the  boundary  conditions  used. 
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Fig.  1 A sketch  of  the  open  and  ducted  rotor  problems 
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Fig.  4 Comparison  of  the  numerically 
obtained  directivity  at  .ft  = 15.0  and  the 
exact  solution  [4]  for  Q = 0.85  with  two 
different  far-field  (PML  and  TAM) 
boundary  conditions  (BD) 


Fig.  5 Comparison  of  the  directivity 
for  different  rotational  speeds 


1 

' (a)  fill 

Fig.  6 (a)  Pressure  lime  history 
sensed  at  ft  = 6.0  and  <t>  = 52°  for 
£2  = 1.15  (b)  Zoom- view  of  the 
box  in  (a):  ducted  rotor  case  with 
M = 0.0. 
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Fig.  10  Comparison  of  noise  radiation  directivity  at  R = 6.0  for  uniform  mean  flow  and 
mean  flow  with  shear  layer 
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Question: 

Do  you  think  this  model  is  applicable  to  swirling  flows  or  only  in  axial  flows? 

Answer: 

Yes.  The  method  itself  is  not  limited  to  axial  flows. 


Discusser’s  Name:  Dr.  Brian  J.  Tester 

Author’s  Name:  Dr.  Yianping  Yin 

Question: 

Is  the  free  shear  layer  of  the  jet  being  modeled  as  well  as  the  shear  layer  on  the  outside  of  the 

duct? 

Answer: 

No.  The  shear  layer  separating  the  jet  and  the  ambient  medium  is  modeled  by  an  analytical, 
parallel  flow  profile.  This  profile  is  only  imposed  on  the  outside  of  the  duct  wall. 


